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Abstract 

Based on the recently proposed Roy-Steiner equations for pion-nucleon {ttN) scattering [T], 
we derive a system of coupled integral equations for the vrvr — > NN and KK NN S'-waves. 
These equations take the form of a two-channel Muskhelishvili-Omnes problem, whose solution in 
the presence of a finite matching point is discussed. We use these results to update the dispersive 
analysis of the scalar form factor of the nucleon fully including KK intermediate states. In partic- 
ular, we determine the correction A^- = a{2M^) — (Jt^n, which is needed for the extraction of the 
pion-nucleon a term from TriV scattering, as a function of pion-nucleon subthreshold parameters 
and the nN coupling constant. 

1 Introduction 

The standard procedure to extract the pion-nucleon a term (TttAt, the nucleon form factor a{t) of the 
scalar current 'm{uu + dd) at vanishing momentum transfer t = (here m = (ttt-u -|- md)/2 denotes the 
average mass of the light quarks), from vrA^ scattering data involves the venerable low-energy theorem 
that relates the Born-term-subtracted isoscalar ttN scattering amplitude at the Cheng-Dashen point 
to a{2Ml) Later on, it was shown that the corrections to this low-energy theorem are very 

small, in particular they are free of chiral logarithms at full one-loop order in chiral perturbation 
theory (ChPT) [HIS]. The extraction of (TttAt itself thus requires knowledge of the difference 

A, = (j{2Ml) - a^N , (1.1) 

which can be determined by means of a dispersive representation of a{t) [6J. The leading contribution 
to the imaginary part originates from vrvr intermediate states, so that, upon neglecting higher terms 
in the spectral function, lTQ.a[t) can be expressed in terms of the scalar pion form factor F^[t) and 
the TTvr NN 5- wave fl{t). Relying on the results of [7J for F^{t) and of [8j for fl{t) led to the 
estimate [6j 

= (15.2 ± 0.4) MeV , (1.2) 

where the error only includes the uncertainty in the parameterization of the vrvr phase available at 
that time. In particular, one should note that the contributions from KK intermediate states in the 
determination of f\{t) and the unitarity relation for a{t) were neglected, while being included in the 
calculation of F^{t). Although the dominant effect may indeed be expected in the pion form factor. 
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such a treatment is strictly speaking inconsistent and leads to an additional uncertainty in (|1.2p that 
is difficult to quantify. Moreover, the result for from [8] corresponds to particular values of vrA^ 

subthreshold parameters and the vrA^ coupling constant. Especially the value 5^/47r = 14.28 used for 
the coupling constant cannot be reconciled with more recent determinations [9Hl2j. 

In this article we perform an updated dispersive analysis of the scalar form factor of the nucleon. 
Based on Roy-Steiner (RS) equations for pion-nucleon scattering [Ij, we determine the tttt — >• NN and 
KK NN S-wmes as solutions of a two-channel Muskhelishvili-Omnes (MO) problem and repeat 
the calculation of the scalar pion and kaon form factors for our input of vrvr and vrvr — )■ KK partial 
waves. Taking everything together, we then analyze the spectral function of the scalar form factor of 
the nucleon fully including the effects from KK intermediate states, and provide an updated value for 
Ao- as a function of vrA^ subthreshold parameters and the ttN coupling constant. 

The paper is organized as follows. In Sect. [2] we consider a generic two-channel MO problem with 
finite matching point and develop a method to construct the corresponding Omnes matrix, which we 
apply to the coupled system of tttt — )• NN and KK — )• NN S-waves in Sect. [3l In Sect. [H we first 
present our results for the scalar pion and kaon form factors, which are then used as input for the 
dispersive analysis of the scalar form factor of the nucleon. We offer our conclusions in Sect. El Several 
technical details of the calculation are relegated to the Appendices. 

2 Two-channel Muskhelishvili— Omnes problem 

We consider the generic coupled-channel integral equation 

. . X 1 /■ , ,T*(t'mt'){(t') 1 f , ,Imf(t') , , 

where bold-faced quantities are two-dimensional vectors in channel space, representing pion and kaon 
intermediate states0 and the imaginary part of f{t) is assumed to be known above the matching point 
tm- In particular, fi{t) and f2{t) can be thought of as the vrvr — t- NN and KK — t- NN S-waves, 
respectively, although the following discussion can be carried out for general angular momentum J. 
The function A(t) contains at most left-hand cuts and is therefore real for t > t-n- = 4M^. The 
unitarity relation is written in the form 

Imf(t) = r*(t)S(t)f(t) , (2.2) 

with T-matrix T{t) and phase-space factor S(t) parameterized a^ 

/r,(i)e2«^^{*)-l \q(t)U^m \ 

= , r?t!w ,We^-(^W-^W)-i ' = diag(arg|^^(t - t.),af fef ^(t - tj,)) , (2.3) 

with pion, kaon, and — for later use — nucleon t-channel momenta 



gt = Jl-M^^ = ^a^, h = Jl-M], = ^af , p, = J\^^ = ^a^, (2.4) 



^Throughout this work, we neglect A-k intermediate states, which are phenomenologically irrelevant at low energies. 
In fact, the effective onset of A-k inelasticities marks the breakdown of the present two-channel model, cf. Sect. [S] 

■^The masses of nucleon, pion, and kaon are denoted by m, M^, and Mk, respectively, and defined by the charged- 
particle masses as given in [13] . 
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and the two-kaon threshold tx = 4M|^. The scattering phases 6{t) and il^{t) are required as input for 
t-K <t < tm, where Watson's theorem [Mj demands ip{t) = 5{t) for t < Ik- Moreover, the modulus of 
the vrvr — >• KK 5-wave g{t) is needed in the full range t-,^ < t < and thus has to be analytically 
continued into the pseudophysical region t-,^ < t < tx. Finally, the inleasticity parameter r/(i) can be 
related to \g{t)\ via 



r/(t) = ^1 - AaJa^{qthYJ\g{t)\'^e{t - Ik) , (2.5) 
and the relation between S- and T-matrix reads 

S{t) = l + 2i Si/2(t)T(t)Si/2(t) . (2.6) 

2.1 Formal solution 

We define the Omnes matrix Q{t) by 

hmnit)=T*{t)mm] ft. <t<tml ^2 7) 

1 Imr2(t) = J 1 otherwise J 

and choose the normalization 17(0) = 1. Writing 

F(t) = f (t) - A{t) = n{t)G{t) , (2.8) 

it follows that for t>t-,^ 

(1 - 2iT* {t)ll{t))n{t+){G{t+) - G(t_)) = 2iT*(t)S(t)A(t) , (2.9) 
where t± = t ^ ie and the physical limit is given by t+ . Using unitarity in the form 

(1 - 2iT*{t)^{t)y^ = 1 + 2iT{t)^{t) , (2.10) 
which in particular holds for t < t/c by virtue of Watson's theorem, we find 

G{t+) - G(t_) = 2in-^{t)T{t)^{t)A{t) , (2.11) 

and thus 

fi(t) f , ,n-\t')T(t')E(t')A(f) n(t) /■ !J-i(OImf(«') 



vr y t' - i vr 7 



dt' 



in 



^ ' TT t' -t TT t' -t ^ ^ 



The problem is such reduced to finding a matrix 0,(t) that fulfills ()2.7p . For 1.,^ < t < tx have 

rj(t+) = (1 + 2ir(t)S(t))l7(t_) . (2.13) 
Taking the determinant on both sides yields, again using Watson's theorem for t < tx, 

detniU) = e^'^^W detO(t_) , (2.14) 
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and thus [15] 

{^m N 
1 I'-JtI^} ■ (2.15) 

Although the determinant allows for an analytic solution in the same way as in the single-channel 
case [16], there is in general no analytic solution for the Omnes matrix itself even for an infinite 
matching point, which therefore has to be calculated numerically, either by an iterative procedure [7] 
or a discretization method, i.e. solving a matrix equation |15lll7j (for a mathematician's point of view 
see [H]). Similarly to the single-channel case, we expect a cusp at ty^, which has to be taken into 
account in the numerical evaluation of the integrals in (j2.12p (see [HIH]). Indeed, for t — )• tm the 
determinant behaves as 

detf7(t)~ |tn,-tr , x = ^^. (2.16) 

vr 

Accordingly, we write for t — ?■ tm from below 

detn{t) = det j7(i„i)e^''^|tm - tr > = ^ij{tm)e'^''^^'^^\tm - ^r'^' , (2.17) 

and from above 

det n{t) = det n{tm)\tra - tr > ^ijii) = ^ij{tm)\tra - > (2-18) 

since i}{t) is real above tm- Here, we have assumed that the (real) functions Ojj(t) are continuous at 
im- The strength of the cusp in each component i^ij{t) of the Omnes matrix is determined by the 
numbers Xij, whose relation to the S'-matrix parameters will be established in the following sections. 
Throughout this paper we will consider the case < Xij < 1, which is relevant for the coupled-channel 
5- wave system of tttt and KK intermediate states. The extension to arbitrary values of Xij can then 
be done along the lines described in [T9|[20]. 



2.2 Dispersive representation of the Omnes matrix 

For < Xij < 1 we may write down a dispersive representation 

where the subtraction constants have been fixed in such a way that 17(0) = 1 and ^{t^) = 0. In 
particular, we can investigate the limit t — )• to obtain information on Xij. Using the asymptotic 
form of the integrals [19] (the dash denotes the principal-value part of the integral) 

*™ oo 

dt' e-s>0,i-s>tm 1 /" dv 

- = — cot nx , 



-t\^ f dt e-S>0,i-S>tm 1 di 

~^ J (t'-t)ltm-t'h W 

tm-e 

dt' e-s>o,f^fm 1 /" 



(2.20) 



vr J {f — t)\tjn — vr J v^{l + v) sinvrx ' 

im 

we obtain for t — )• from belo'vJl 

smvrxii 2i sinvr3;2i 



^We consider the case J = for simplicity. The general case can always be recovered by introducing the correct 
phase-space factors according to cr^ — > o"??*'' and — af'kt''. 
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i'i2 V — t e — iii2 tm — i e ^ — l-5i22 rm — S e ^ \g\at e 

sin7rxi2 2,1 sin7rx22 

sinvrxii 

„i7rX2i 1 _ „p-2j(i/)-5) 
+ f^2l|tm - tr-e^^- ^ , 

sm 7rx2i 2« 

%7TX X2 

sin7rxi2 

e*'^^22 1 _ ^( 

sin 7rx22 2i 



i7rx22 1 _ „p-2i(i/)-5) 

+ ^^22|tm - tr^^e^^^^-^ , (2.21) 



and for t ^ t^a from above 

piSu 1 _ „p— 2j5 _ p«'52i 

Oiiit^ - tr- = oiiitn. - tr^-. — + n2iK - tr^-. — lowte-'^ , 

sinvrxii Zi sin7rx2i 
piSi2 1 _ ^p—2iS _ pi522 
^12\tm - = n^ltm - tr'-. + f^22|tm " t^^' -. l^kf^e"*'^ , 

Sin7r3;i2 2i Sin7rx22 

piSii _ pi&2\ 1 _ ^p-2i{ip-S) 

^2i\tm - = ^uK - r'"^ bk^e"^^ + n^iK - , 

sinvrxii sin7rx2i 2i 

piSl2 _ piS22 1 _ j,p-2i{ip-5) 
^22\tm-tr^ =nu\tra-tn-. 1 e"^^ + 1^22 1 " C ^ , (2.22) 

Sin7r2;i2 sin7rx22 2^ 

where we have suppressed the evaluation at wherever possible. If we assume that g(tm) 7^ and 
^ij{tm) 7^ (which can always be achieved by choosing the matching point appropriately), we can 
conclude from the first line of (j2.2ip that X21 > xu, since otherwise g or CI21 would have to vanish 
at tra- Conversely, the third line requires 2:21 < xn by the same argument, and hence xn = X2i- 
Similarly, we find 2:12 = 2:22 • Moreover, as the determinant behaves according to (|2.16p . we can 
conclude that xu + 2;i2 = x, again provided that detil(tm) 7^ 0. Dividing the first line of (12.2ip by 
the first fine of ^22i), we find 

]=e*™2i[i_^ \ (2.23) 

sinvrxii 2i J \ sinvrxn 2i J 

which for 2:21 = xu reduces to 

and thus nxu = 5u up to integer multiples of 2-ir. Arguing analogously for X12, these results can be 
summarized as 

S 

Xu = X21 , 2:12 = 2:22 , Xij = — , xii + 2:12 = 2; . (2.25) 



By virtue of K2E\t . K2T\i and ^^2% take the form 



smTTXii 2i sinirxiil^i tm 11 

K sinTTXii I "I im sinTTXii 2i / \ 

^ 2i ^ sin7r.T,ol5'Fi„e \ ^12 



21, 



sin7rxi2 2i sin7rxi2 1^1 *m 

e''^^l2 I I TT p-itl> e''^^l2 l-r?e-^'(^-^) 



^ \q\ai-' e — tt- i/ \f^22 

sm7ra;i2i^i sm7rxi2 2i / \ 



. (2.26) 
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To ensure the existence of non-trivial solutions the determinants of the coefficient matrices must 
vanish. This leads to 



cos7r(2a;ii— x) — r/cos7r(2y — x) = , cos7r(2xi2 — a;) — r/cos7r(2y — x) = 



V 



Sit. 



IT 



(2.27) 



These conditions are invariant under xu — t- x — xu, i.e. there is an ambiguity between xn and xi2. 
However, demanding that xn coincide with x in the single-channel limit yields 



Xii 



-|x H — arccos(7y cos vr(2y — x))| , xi2 = -|x 
2 I vr J 2 1. 



1 



arccos 



TT 



(r/cos7r(2y-x))| . (2.28) 



Finally, (l2:26]l together with (|2:28|) and 



sm TTXii = sm ■ 



TTX 1 + Z 



+ COS 

2 2 



TTX 1 — Z 



sm 7rxi2 = sm 



vrx 1 + z 



TTX 1 — Z 



COS ■ 



COS TTXll = COS ■ 



vrx 1 + z 



TTX 1 — Z 



2 V 2 

z = ??cos7r(2y — x) , 

can be used to derive constraints on Q,-,. We find 



sm ■ 



TTX 1 + Z TTX I — Z 

cos7rxi2 = COS— W^— -hsm— W 



(2.29) 



21 



N 



^11 'i\9Wt 



^22 



N 



N = \J\ — T]^ cos^ '?r(2y — x) — r/ sin vr(2y — x) . (2.30) 



In the single-channel case one can show that, using the integrals (j2.20p . the solution for /(t) is 
automatically continuous at t^^ [19j. In fact, the same statement holds true also in the two-channel 
case. The relations (j2.30p are essential in the proof, as demonstrated in Appendix [Al 



2.3 Construction of the Omnes matrix 
2.3.1 Infinite matching point 

Our construction of the two-channel Omnes matrix with finite matching point will heavily rely on the 
solution for its infinite-matching-point analog r2°°(t), whose defining property can be stated as 



Iml7°°(t) = r*(t)S(t)S7°^(t) 
Iml7°^(t) = 



for 



1 otherwise j 



For its calculation we follow [T5] and discretize the unsubtracted dispersion relation 



(2.31) 



' vr J t> -t 



(2.32) 



on a set of Gaufi-Legendre integration points. Note that an unsubtracted dispersion relation converges 
provided that the phase-shift at infinity is positive. In the one-channel case this can be directly deduced 
from the explicit solution 



{LXJ 

t-TT 



(2.33) 
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which behaves as 

(2.34) 

for large t. The unitarity condition (|2.3ip can be rewritten as 

' rism{25-ip)+sm'tp 2|g|o-f gft-ty ) 

imnr = I ''<=-(2t,ti^-«'^ "Z^^^J:^!. I Renr , ^ e {1,2} , (2.35) 



with 



2\g\crf ri sm{25—ip) —sin ip 

T] cos(25— i/))+cos 1/1 cos(2(5— i/i)+cos -ip 



^oo _ I 1 000 _ ( ^12 



= , n~ = , (2.36) 



which below the two-kaon threshold reduces to 

tan 6 

cos 5 



Imnr={ ijwi Re^r- (2.37) 



The details of the numerical solution of the corresponding integral equation for Re are described 

in mm. 



2.3.2 Finite matching point, single-channel case 

An Omnes function with a finite matching point does not allow for an unsubtracted dispersion relation, 
since the solution (j2.15p tends to a constant for t — )• oo, and consequently one picks up contributions 
at infinity. Moreover, the cusp at the matching point renders both the discretization method and 
an iterative procedure involving subtracted dispersion relations inappropriate, as neither is able to 
accurately reproduce the analytic behavior around t^a- For this reason, the aim of this section is to 
establish a method that relies on the known solution in the infinite-matching-point scenario. 
We first observe that the function 



x{t) 

m 



x(t) 

^~nvx(t)e(t-t,-n) (2.38) 

has the correct properties to cancel an imaginary part above tm- Indeed, the function 

J7(t) = !^°°(t)e(t) , (2.39) 

with f]°°(t) from (ITHHI) . fulfills 

|lmJ](t) = J]°°(t)|e(t)|(e*('^^W-^W)sin(^(t) -sin7rx(t))j \ t>t^ ]' ^' ' 

and with the choice 

for <^ - (2.41) 





the defining property ()2.7p holds. Since we know the analytic solution, we can study the properties of 
this construction in more detail 



t'TT tu 



7 



The first term coincides with the expected result for the single-channel case, while the second factor 
is new. It is indeed real, and thus preserves all defining properties, in particular the normalization 
r2(0) = 1. This example shows that we would exactly recover the result (j2.15p if we did not know the 
solution for a finite matching point, constructed it according to (|2.39p . and chose 5{t) = 5{tm) for t 
above the matching point. Obviously, ()2.42p implies that "the" Omnes function is not unique. The 
derivation of (I2.12p . however, only relies on the defining properties and is therefore independent of 
such modifications. 



2.3.3 Finite matching point, two-channel case 

To generalize the preceding considerations to the two-channel case we define a matrix 



^^At) = ) (2-43) 

with real functions iij{t). In view of the results of the previous section we take phases and inelasticities 
constant above and thus can ignore the t-dependence of Xij and which in the following will 
always be understood to be evaluated at tm. It is straightforward to show that 



m = a(t)(!!~)'^ ((){{() I 1 1 > !,„ I ' 



t t^a 



with infinite-matching-point solution and a real matrix a{t), fulfills (j2.7p provided that 

Im^^(t) + ^^(i)r(t)S(i) =0 for t > , 

[a{t),T*{t)T.{t)] =0 fort<tm. (2.45) 

Imposing xu = X21 = x — X12 = x — X22, the first condition corresponds to 



sinTrxii 1i sin-Trxn 1^1 








sin7ra;i2 2i sin7rxi2 1^1 *m 1 ( "^^^ \ _„ (0 Ae(\ 

sin 71X12 I "I *m sin 71X12 2j 

Non-trivial solutions of (|2.46p again exist for xu and xi2 given by (j2.28p . while the components of 
^(t) are related by 

6i N 62 N 



For definiteness, we take 



which ensures that 



61 m^y^ 62 2\g\al 



K 



(2.47) 



61 = 1 , 62 = ( 1 + 2 ) ' (2.48) 



det^(t) = (^^) , (2.49) 



I'm 



and thus, by the results of the previous section, would preserve the form (j2.15p of the determinant of 
the Omnes function if deta(t) = 1. The condition (j2.45p on a{t) requires 

aj r] sin(2(5 — ip) 

021 =012—^! 022=011-012 --—77 , (2.50) 
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while Oil and ai2 can be chosen freely. 

Finally, i}{t) should fulfill the normalization Q{0) = 1, which can be achieved by modifying the 
normalization in the calculation of r2°°(t) appropriately. Assuming a(0) = 1, the corresponding 
condition 

e = - l> - =1 2.51 



leads to 



n-io) = ( 1 + 1 I _.r N "p:^ I • (2-52) 

The above construction (j2.44p ensures that ^ij{t) has the expected behavior for t — )• tm, that the 
factors of (1 — t/tm)^'^ factorize in (I2.46P (which have therefore already been canceled there), and that 
detQ(t) is continuous at tm. However, in general, Clij{t) itself will not be continuous at and the 
condition (j2.30p will be violated. 

In order to remove these shortcomings we make use of the freedom in choosing a{t). In fact, a 
particular choice of a{t) can enforce continuity of either fii(t) or f^2(i)) but not of both simultaneously. 
This impediment can be circumvented by noting that (I2.44p may be regarded as separate equations 
for fii(t) and Cl2{t). We can thus derive an 0!"^\t) from a construction with an ai{t) tailored for this 
component (discarding Q!^\t) in this case), Q!^\t) from a different a2{t) (discarding ^^^\t)), and 
finally join these two vectors into the final Omnes matrix Q.{t) = {f2^^'*(t), il^\t)^. 

Below the two-kaon threshold we take ai{t) = a2(t) = 1, while for t >tK 



aiit) 



1 f(t) 



1 -w/(0 

\ N> 

with 



N{t) = v"! -r?(i)2cos2(2(5(t) - ^(t)) - r?(t) sin(25(t) - tpit)) , (2.54) 

cf. (|2.30p . and a function /(t) fulfilling fitx) = 0, /(tm) = proves adequate. This construction 
makes sure that ^ij{t) is continuous at and that the relations ()2.30p hold. 

As we have seen in the previous section, the Omnes function is not unique, and therefore there 
is a priori no reason why the determinant of the resulting 0(t) should match the single-channel 
expectation: we can always multiply ^i{t) with a real function g{t) with g[{)) = 1 without vitiating 
the above construction. Since the unitarity condition in its form analogous to ()2.35p alone implies 
that 

det Q. = — (Re OnRe O22 - Re Jli2Re O21) , (2.55) 

•q cos{2o — ip) + cos yj 

and thus ensures the correct phase, this could be used to make the modulus of the determinant coincide 
with (|2.15p . We simply take 

This choice of the exponent ensures that f{t) decreases rapidly below tm, but not so fast as to 
cause numerical problems in the matching-point region. In this way, the single-channel result for the 
determinant is accurately reproduced except for the energy region close to t^. 
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3 Application to pion— nucleon and kaon— nucleon scattering 



The vrvr — )• A^A^ and KK — )• A^A^ 5"- waves /^(t) and h\{t), respectively, fulfill the unitarity relations 
(cf. [1] and Appendix iBj) 

lm/^0(t) =cTf (rO(t))*/.°(t)^(t-tx) +^cTr(5o°(t))7^(t)^(t-g , (3-1) 

where t^, g^, and rg denote the vrvr, vrvr — )• KK, and KK 5"- waves with t-channel isospin It = 0. Since 
by virtue of unitarity in the vrvr / KK system we may identify 



4it) = ' , = b(t)|e^^(*) , = ^^^^^^^ , (3.2) 

the unitarity relation p.ip reduces to 



imf(t) = T*{t)mm , m = ( ' (^.s) 

and with RS equations providing a dispersion relation of the form 

oo 

m = Mt) + (a + b*)« - + / dt' ,,,,*^^;;;,,,) . (3.4) 

where t]\f = 4m?, the discussion in Sect. [2] applies. We consider here only the twice-subtracted version 
of the RS equations for vrA^ scattering derived in [1]. The corresponding inhomogeneity for the vrA^ 
system reads 

A^it) = Alit) = Nlit) + Alit) , 

1 

Alit) = - J dTy'^{Goi|'""^t,W^Olm/,+ (M^O + Goi|'""'^t,-W^Olm/(|+i)_(M^')} • (3-5) 



(=0 



Here, N^{t) denotes the nucleon pole term, fl^{W) the s-channel partial waves with total angular 
momentum j = |/ it 1/2| and isospin index 7 = it, and W+ = m + M-,^. Moreover, we have already 
neglected the contributions from higher t-channel partial waves, which were shown to be insignificant 
in [1]. In fact, the corresponding equation for A2{t) is very similar, as long as we restrict ourselves to 
combinations of isospin and angular momentum that couple to vrvr (cf. Appendix [B]): the nucleon pole 
terms need to be replaced by the hyperon-pole contributions (|B.10p . by Mk-, and the Clebsch- 
Gordan coefficients in the relation between the s-channel amplitudes in I = ± and = 0, 1 bases 
are different from vrA'^ (where Is = 1/2, 3/2). In particular, the kernel functions GQi{t,W) require no 
further modification besides M^- — )■ Mk- Finally, the subtraction constants ai and bi can be related 
to standard vrA^ subthreshold parameters 

ai = -— ( - + 4n + ^nn— V ^1 = -— Mm"— 1 (3-6) 
levr I m °° °° 3 i 16vr I °^ 12 i ^ ' 
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(see [T] for precise definitions), and similarly for 02, 621 and KN subthreshold parameters. 

Starting from the dispersion relation (I3.4p . the solution for f(t) can be derived along the lines that 
led to (j2.12p . In addition, we may use the spectral representation of the inverse of the Omnes matrix 



to perform the integrals involving a and b explicitly (J7 denotes the derivative with respect to t). In 
this way, we arrive at 

f(t) = A(t) + {t- tNW){l -tn{0))a + t{t - tN)^{t)h 



TT ^'J t'^t' - tN){t' - 1) vr 'V t'2(t'-tiv)(i'-0 ' 

To obtain the corresponding numerical results for /^(t) and h^{t) we now collect the various 
input needed for the explicit evaluation of (j3.8p . We take the vrvr phase and inelasticity from the 
extended Roy-equation analysis of [21''22] and the tttt — >■ KK partial wave from [19], where, in the 
pseudophysical region t-,,^ < t < tx, the modulus \g{t)\ was determined as the solution of RS equations 
for ttK scattering, while above the two-kaon threshold phase-shift solutions |231I24| were usedQ The 
two-channel approximation in terms of vrvr and KK intermediate states should work well at lower 
energies where the /o(980) resonance dominates, but will break down once inelasticities arising from 4vr 
intermediate states become important, which are expected to set in around ^/Iq = 1.3 GeV. Evidently, 
this still leaves a large part of the vrvr — )• NN pseudophysical region t-j^ < t < tjy where intermediate 
states besides those considered here may contribute significantly to the spectral function. However, 
two subtractions should already appreciably suppress the sensitivity to this high-energy regime, which 
we will examine in more detail below. 

In view of the remaining uncertainties in the pseudophysical region we refrain from including 
information on f^{t) even above tN from a phase-shift solution [25], but just put Imf = above 
tra- Due to continuity of the MO solution at tm this implies that the solutions for and h^{t) 

will vanish at tm as well (unless the phases happen to pass an integer multiple of vr at tm)- We take 
advantage of the fact that for kinematic reasons (t) and (t) have a zero at the physical threshold 

and choose the matching point accordingly as tm = ^n, which should allow for a reasonably smooth 
matching. 

Next, to evaluate the inhomogeneities Aj(t) we need the imaginary parts for the s-channel partial 
waves of vrA^ and KN scattering. As discussed in [l], we rely on the GWU solution [26l[27] for 
W < Wa = 2.5 GeV in the vrA^ case, summing up all partial waves up to /max = 4, and on the Regge 
model [28] for W > Ws.- Similarly, A2 is evaluated based on [27^129]. integrating the partial waves 
with Z < 4 up to 2 GeV. 

In the remainder of the paper we will consider three variants of the input described above. First, 
we take 5 and ^ to be constant above to ("RSI"), second, we guide 6 and ip smoothly to 2vr above to 
by means of pTS] 

5{t) = 2n + (5(to) - 2vr)/f A) , />) = , (3.9) 



keeping the phase constant above tm ("RS2"), and third, we modify RSI in such a way that A2 = 
("RS3"). In all three cases, \g{t)\ is led smoothly to zero above tm by a prescription similar to ()3.9p . 



We are indebted to Bachir Moussallam for providing a version of the solution for g{t) consistent with the tttt phase 
shift of [2n[22] . 
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Figure 1: Results for real and imaginary part of the components of the Omnes matrix for the input 
phases RSI and RS2 as described in the main text. 



The choice of these variants is motivated as follows. As indicated above, the model for the ttti/KK 
^-matrix is only meaningful roughly up to Iq, and ideally our results should be insensitive to variations 
of this input above Iq, the simplest choice of course being to keep the phases constant. However, to 
ensure the correct asymptotic behavior for the scalar meson form factors, the phase ip must tend to 
an asymptotic value of 27r [15], which, phenomenologically, also suggests to guide 5 to 27r. Thus, RSI 
and RS2 are convenient choices to assess the sensitivity to the high-energy input for the phase shifts. 
The results for the Omnes matrix corresponding to these two scenarios are depicted in Fig. [TJ In 
addition, we compare the results for detr2(t)/(l — t/t^Y to the analytically known results in Fig. [2l 
As expected, the only deviations occur in the proximity of tnu where the modifications originating 
from f{t) according to (I2.56P set in. Note that the difference between RSI and RS2 appears slightly 
exaggerated here, since x differing in both cases leads to a different factor being divided out. 
Finally, we use the reference point [8] 

1^ = 14.28, 4o = -l•46M-^ d+=1.14M-3, 6+=-3.54M-3, (3.10) 

although we will display the dependence of Ag- on each of these parameters explicitly in the end. In 
contrast, we are not aware of reliable input for the KN subthreshold parameters, and simply put 
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Figure 2: Comparison between numerical and analytical result for the determinant of the Omnes 
matrix. 

02 = 62 = 0. In fact, this approximation gives reason to investigate RS3, since the results from the 
vrA^ sector show that the contributions from the corresponding parameters will certainly not be larger 
than the sum of KN Born terms (conventions for masses and couplings of the hyperons are given in 
Appendix [B|) and s-channel integrals. The difference between RSI and RS3 thus serves as an estimate 
of the uncertainty induced by neglecting the KN subthreshold parameters. 

The results for /°(t) and h\{t) are shown in Fig. [3l In the case of /+(t) the agreement between 
the different parameter sets is very good up to 1 GeV. It is striking that the difference between RSI 
and RS2 is much larger than between RSI and RS3, which indicates that the results are much more 
sensitive to the choice of the phases beyond to than to the details of the KN amplitude. The real part 
of /i4.(t) exhibits two distinct divergences below 1 GeV that correspond to the pole-term contributions 
from the A (large peak at 0.82 GeV) and S (small peak at 0.73 GeV) hyperon, strictly analogous to 
the nucleon pole in /+(t) that emerges slightly below t^^. These poles, which disappear in RS3 by 
construction, do not pose a problem in practice, since /^^(i) only contributes to the spectral function 
of the scalar form factor of the nucleon above Ik- Apart from these poles below the two-kaon threshold 
the conclusion is very similar to f^{t): the uncertainty in the phase shifts outweighs the uncertainty 
in the KN input. 

4 Scalar form factors 

4.1 Scalar pion and kaon form factors 

We define the scalar pion and kaon form factors as 

F^{t) = {7r{p')\m{uu + dd)\7r{p)) , F|(t) = {K{p')\m{uu + dd)\K{p)) , t = {p' - pf . (4.1) 
In the two-channel approximation they fulfill the unitarity relation [7] 

ImF^(t) = T*{t)m^'{t) , F^(t) = f^ll^)) ' (4-2) 

and thus 

¥^{t) = an?° + /3n^ , (4.3) 
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Figure 3: Results for real and imaginary part of and The solid, dashed, and dot-dashed 

lines refer to the input RSI, RS2, and RS3 as described in the main text. For the black crosses 

indicate the results of [8]. 



with the infinite-matching-point Omnes solutions fl^ defined in ()2.36p . The phases 6 and ip are 
guided smoothly to their assumed asymptotic value of 27r according to ()3.9p . Using the normalization 
r2°°(0) = 1 of the Omnes matrix to pin down the coefficients a, /3, we find 

Fkt) = ^F^mTiit) + FimT2{t) . (4.4) 

The form factors at t = can be determined via the Feynman-Hellmann theorem |30y31j 

^'^^^ = ' = ^^^^ ' ^^-^^ 

from the quark-mass dependence of the meson masses. For the pion form factor the result at 0{p'^) 
in the chiral expansion reads [32j 



F^iO) = Ml - ^^{h - 1) = (0.984 ± 0.006)M2 , (4.6) 
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OA Ml 


Af2/2 


0.6 Ml 


(r2)| [W] 


0.575 
0.835 


0.584 
0.710 


0.592 
0.626 



Table 1: Dependence of the scalar pion and kaon radii on FJ-(O). 



where we have used Z3 = 3.2 it 0.8 [33]. The leading-order result for the kaon form factor 

]\/f2 

F|(0) = ^ (4.7) 
is subject to potentially large SU{3) corrections, which in the isospin limit amount to |34j 
o M2 f Ml Ml Ml r Ml \ m2 M2 

^'-^^^(2Li-m^'-^(2Li-L-A, (4.8) 

where is the mass of the rj, the pion decay constant, and fi the renormalization scale. Varying 
the low-energy constants in the range |33j 



2L^-L^ = (-0.35...+0.1)-10-^ , 2L^ -L^ = (0... + 0.2) • 10-3 , (4.9) 

corresponds to F^{0) = {OA. . . 0.6) M^. In the following, we will restrict F^{0) to lie within these 
boundaries, while adopting F^{0) = Ml/2 as our central solution. 

The corresponding results for the form factors are depicted in Fig. [H The phase of coincides 
with S below tx, as required by unitarity, cf. (|2.37p and (|4.4p . Above tx, the behavior of the phase 
actually depends on the assumption for F^{0), it either largely follows S or abruptly drops by tt. 
The appearance of the first scenario is surprising in view of the results of [35] (see also [36j), where 
only the second behavior is mentioned (the assumption for F^{0) agrees with our central solution). 
The reason for this discrepancy can be understood as follows. The 5-matrix used in [35] involves a 
phase 6 that fulfills Sitx) < vr. In this case, the phase 6t of the full vrvr partial wave ig itself displays 
the characteristic drop above tx, reflecting the fact that the phase arrives at the two- kaon threshold 
immediately before completing a full circle in the Argand diagram. However, in recent years, it seems 
to have become consensus that 6{tK) > vr is more likely [211122^1371138] . which implies that 6t by no 
means exhibits a sharp drop above tx- We conclude that the behavior of the phase cannot simply 
be deduced from the phase of tg, it crucially depends on the relative strength F^{0)/F^{0) in the 
superposition of the two terms involving different components of the Omnes matrix as given in (j4.4p , 
thus attesting to the inherent two-channel nature of the problemU Prom this point of view, it is not 
surprising that the phase may behave differently if F^{0) is varied. In contrast, the phase of F^{t) 
roughly follows the shape of ip for all three solutions. Note that in this case (|2.37p does not impose 
any additional constraints on the phase below tx- 



^One immediate consequence is that an effective single-channel Omnes description of in terms of the phase of 
- only be applic, 
term —nO{t — Ik)- 



will only be applicable for certain ranges in F^ (0) /F^ (0) , unless the phase is supplemented by hand with an additional 
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Figure 4: Modulus (left) and phase (right) of the scalar pion and kaon form factors. The solid, 
dashed, and dot-dashed lines refer to F|(0) = Ml/2, OAM^, and 0.6 M^. The phases of and F| 
are compared to 5 and ip, respectively, as indicated by the dotted lines. 



Finally, we can express the scalar radii in terms of the form factors at t = and the derivative of 
the Omnes matrix 

(r2)f = 6|r2??(0) + -^||||iir2(0) 

(^')i = 6{^^^^(0) + f^^2(0)} . (4.10) 
In this way, the derivative of the Omnes matrix, e.g. for our central solution 

r2-(0) = (^-^^ ^■^^] GeV-2 , (4.11) 
\1.26 0.89 J 

leads to the results for the scalar radii summarized in Tabled) Our results for the scalar pion radius are 
in good agreement with (t'^)^ = (0.61 ± 0.04) fm^ from [39] and the range (r^)^ = (0.583 . . . 0.653) fm^ 
found in [15] . Albeit attached with a fairly large uncertainty, the values for {r"^)^ lie systematically 
higher than its ChPT expectation (r^)^ ~ 0.3 fm^ [20] (for a detailed comparison of the dispersive and 
the ChPT result as well as the role of 0{p^) corrections see [S]). In both approaches the uncertainties 
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(FiiO), A) 


{OA Ml, 1.3 GeV) 


{Ml/2, 1.3 GeV) 


{0.6 Ml, 1.3 GeV) 


{Ml/2, 2m) 


RSI 


13.81 


13.86 


13.91 


13.65 


RS2 


13.75 


13.89 


14.04 


13.56 


RS3 


13.92 


14.00 


14.09 


13.80 



Table 2: Results for A^- in MeV for various combinations of F^{0), integral cutoff A, and input. 

are substantial, either due to the large sensitivity to the specific input in the dispersive calculation 
or due to insufficient knowledge of low-energy constants and higher-order corrections. As the precise 
value of the scalar kaon radius is irrelevant for the present study, we do not consider this issue any 
further. 

4.2 Scalar form factor of the nucleon 

The scalar form factor of the nucleon is defined as 

= ^{N{p')\m{uu + dd)\N{p)) , t = {p'- p)2 . (4.12) 
Zm 

It fulfills the once-subtracted dispersion relation 

oo 

= + dt'^, (4.13) 

where the pion-nucleon a term cTttAt = cr(0) acts as subtraction constant. In this way, evaluation at 
t = 2MI gives access to Ao- (as defined in (jl.ip ). provided that the imaginary part is sufficiently well 
constrained to perform the dispersive integral. Generalizing the result quoted in [6] by including KK 
intermediate states, the spectral function becomes 

lma{t) = --l^|^g,(Fi^(t)) VO(t) e(t - g + h{F§{t)yhl{t)e{t - t^) | . (4.14) 

The corresponding results using the input RSI, RS2, and RS3 as discussed in Sect. [3] as well as 
our central solution for the scalar pion and kaon form factors are depicted in Fig. [5l We also show 
a variant of the spectral function where the second term in the above unitarity relation (|4.14p due 
to KK intermediate states is neglected. While the impact on the real part is moderate, we see that 
the spectral function develops an imaginary part starting at tx- In contrast, our full solution stays 
real as long as the input for the phases is treated in the same way in the calculation of the meson- 
nucleon partial waves and the scalar meson form factors. For this reason, the results for RSI and RS3 
become complex around to, while RS2 is real in the full energy range (apart from some numerical noise 
at the two-kaon threshold). These findings emphasize the importance of treating inelastic channels 
consistently in all contributions to the unitarity relation, in particular the necessity to explicitly include 
the intermediate states that are responsible for the inelasticities. 

The dependence of Ag- on the cutoff A of the dispersive integral is shown in Fig. [6j We see that 
the dispersion relation converges quickly and the results hardly change above 1 GeV. We quote the 
outcome for A = ^/to and input RS2 as our final result, estimating the uncertainty by the variation 
induced by changing the cutoff to 2m and varying the input set or the assumption for F^{0) (cf. 
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Figure 5: Real and imaginary part of \ra.a{i) for our full solution (top) and with set to zero 

(bottom). The solid, dashed, and dot-dashed lines refer to the input RSI, RS2, and RS3 as described 
in Sect. [3l 

Table [2]). Moreover, /+(t) and /i'^(t) depend linearly on the vrA^ parameters (j3.10p . so that the 
corresponding corrections for changing these parameters can be determined straightforwardly. Putting 
everything together, we find 

= (13.9 ± 0.3) MeV 

+ Zi (^1^ - 14.28^ + Z2 (4o + 1-46) + ^3 (4i Ml - I.m) + Z4 (&+ ^1 + 3.54) , 

Zi = 0.36 MeV, ^2 = 0.57 MeV, Z3 = 12.0 MeV , ^4 = -0.81 MeV. (4.15) 

These results are in reasonable agreement with [6], and remarkably close to the O(p^) ChPT analysis 
of [12] (for earlier work in ChPT on see [iHtiiS] ) 

= 14.0 MeV + 2M^e2 , (4.16) 

where 62 is an 0{p'^) low-energy constant^ The potentially largest correction in (j4.15p originates from 
d^i, e.g. taking d^^ = 1.27 from [37] increases Ag- by 1.6MeV. In contrast, adjusting the coupling 
constant to g'^/^ir = 13.7 [111 112] only leads to a correction of —0.2 MeV. Indeed, this result is not 
surprising, as it is that controls the slope of the scalar form factor of the nucleon. 

®In fact, 62 contributes to the chiral expansion of djj. 
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A / GeV 

Figure 6: Ag- as a function of the integral cutoff A. Line code as in Fig. [5l 

5 Conclusion 

We have presented a method to construct two-channel Omnes functions for a finite-matching-point 
setup, and applied this formalism to solve a coupled-channel integral equation for the tttt — )• NN and 
KK —7- NN S-waves, which can be derived from Roy-Steiner equations for vrA^ and KN scattering. We 
have also repeated the conventional two-channel Muskhelishvili-Omnes calculation for the scalar pion 
and kaon form factors, which together with the meson-nucleon partial waves determine the spectral 
function of the scalar form factor of the nucleon fully including the effects of KK intermediate states. 
Based on these results, we have updated the dispersive analysis of the correction A^- that relates the 
pion-nucleon a term to the vrA^ amplitude at the Cheng-Dashen point. Throughout the calculation we 
have investigated the sensitivity to various pieces of the input in detail, in particular, we give our final 
result for as a function of ttN subthreshold parameters and the vrA^ coupling constant. Our result 
essentially confirms the result of [6], however, using modern phase-shift input and with uncertainties 
due to KK effects much better under control. 

Besides its implications for the scalar form factor of the nucleon, the present calculation of the 
vrvr — )• NN S-wave in a full two-channel treatment completes the solution of the t-channel part of the 
Roy-Steiner equations derived in [1], and is also a key ingredient to pin down the vrA^ amplitude at 
the Cheng-Dashen point [HEKH]. Moreover, the formalism developed here should prove valuable for 
other systems requiring finite-matching-point two-channel Omnes solutions as well, e.g. for including 
KK effects into the Roy-Steiner analysis of 77 — > tttt [20] . 
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A Continuity at the matching point 

To prove the continuity of f(t) at tm, we rewrite (j2.12p in terms of a principal- value integral as 



f(t) = (1 + ir(t)s(t))A(t) + — dt' — y ^ y ' y ^ v ; _^ _w / 

TT J t — t J 



t' -t 



■ (A.l) 



We only consider the case t — )• from below (continuity from above can be proven in a similar 
way). In this limit, the whole mass of the integral is concentrated at tm- Using (12.17p . (I2.18p . (12.20p . 
and (|2:25]) . we find 



f (V) = (1 + i r(V)E(t„,)) A(U) + 



h hi 



detO(tni) Win ^iv 



Tit^)^it^)A{tr 



1 



h hi 



detil(tm) Win -^iv. 



Imf(tm) 



(A.2) 



with 
h 
hll 



011^22 cot vrxii — rii2^^2i cot 7rXl2 ) e 



ill) 



^^21^^22 ( cot vrxii — cot 7rxi2 ) e 



hi = r2iirJi2 ( cot vrxii — cot 7rxi2 I e 



^11^^22 cot 7rXl2 — $712^^21 COt TTXn ) 6 



/i = — r2iir222- 



sm vrxii 



+ r2i2f^2i" 



sm 7rxi2 



/ll = Oiiil 



12 



oi'5i2 



smvrxii sm7rxi2 



hll 



-r22lf^22 



smvrxii sm7rxi2 



IlV = —^^11^^22- 



o«<5ii 



sm 7rxi2 



+ 17x2^^21" 



smvrxii 



. (A.3) 



For the continuity condition (|A.2p to be fulfilled, in particular all terms depending on A(tni) must 
cancel amongst themselves. Putting the coefficients of A,(tjn) in each component fi to zero yields four 
constraints on h^hv, which can be inverted to obtain 



h = idetCl e* 



hv = i det Q 



B 



/ii = ^/iii = 4detfie 
(J 



2iivx 



K 



(A.4) 



^ + g2i7rx^* ' 

where 

A = r^e^^^^'^-y) - 1 , B = r^e^^^^ - 1 . (A.5) 

Similar considerations apply to h~hv- By means of Imf = r*Sf, the remaining pieces of ()A.2p can 
be expressed as linear combinations of /i and /2, which altogether again amounts to four constraints. 
The solutions are 



h 



-2i det e 



2inx 



A + e^^'^^A* 



hy 



-2i det e 



2j7rx 



(T- 



0-;' 



K 



-hll = hie 



(A.6) 

Based on (j2.29p it is straightforward to show that ()A.4p and ()A.6P hold as long as the following 
relations are fulfilled 



\/\ — (f^iif^22 + f^i2f^2i) = Ti sin 7r(2y — x) det VL , 

l^lfif" detJ^ , 



- z20iil7i2 = -blo-/^ detJ^ , Vl - z'^^21^22 = l^ki''^ detl^ . (A.7) 

In this way, the properties (|2.30p of the Omnes matrix finally ensure continuity at the matching point. 
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B Conventions for kaon— nucleon scattering 

In general, we follow the notation and conventions for vrA^ scattering of [Ij. In this appendix, we briefly 
collect the additional formalism necessary to arrive at a closed system for vrA^ and KN scattering. 
We define the kaon states 



1 1 

2'2/ ' 



1 1 
2'~2 



1^" 



1 1 

2' 2 



with crossing properties 



C\K+) = -\K-) , C\K^) 
and choose the isospin decomposition of the amplitude as 



1 1 

2'~2 



(B.l) 
(B.2) 



T = T+ -tn-tkT- =T+ -2[h{h + l)--)T~ 



(B.3) 



which leads to 



rj-ilu=0 3T~ 

rplt = _ 27''^ 



rpls = l _ rp+ rp- 



rplu = l 
rplt = l 



2T- . 



(B.4) 



In these equations Ig, It, and denote s-, t-, and u-channel isospin, respectively, and tn and tk are 
the Pauli matrices associated with the nucleon and kaon isospin operators. In particular, we have 



(B.5) 



In complete analogy to ttN scattering, the amplitude for the process K(q) + N(j)) — )• K(q') + N{p') 
can be decomposed as 



u{p) , /G{+,-} 



(B.6) 



The Born-term contributions due to hyperon pole diagrams arq3 



iK+p 



E2 my — m 
9kny 



y=A,s 



u — m 



2 ' 
Y 



y=A,s 



s — m 



2 ' 

y 



B 



K+p 



E 

y=A,E 



9kny 



u — m 



2 ' 

y 



qK-p^ 



9kny 



y=A,s 



s — m 



2 ' 

y 



(B.7) 



and hence 



— ^ 9kny 



my — rn 



y=A,E 



1 1 

± 



s — u — my 



y=A,s 



9kny 



=F 



s — my u — my 



^We use niA = 1.116 GeV, ms = 1.193 GeV [13], and Qkna/'^''^ = 15.55, Qknt./'^''^ = 0.576 (48j for the masses and 
couplings of the hyperons, respectively. 
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The KK — )• A^A^ partial waves can be obtained from the invariant ampUtudes by means of the 
projection formula (with the t-channel scattering angle zt = cos 9^^) 



/i:J^*=0'i(t) = -^ I dztPj{zt] 



■iTT 



A^{t,zt) 



m 



{pthy- 



-ztB^{t,zt) 



-1 

The main difference to the vrA^ f-channel partial-wave projection [39] originates from the fact that 
due to the lack of Bose symmetry in the KK system a partial wave with given angular momentum 
J couples to both It = and = 1 (corresponding to + and — on the right-hand side of ()B.9P ). In 
the following, we are only interested in the combination where even/odd J corresponds to It = 0, 1, 
respectively, since only these partial waves can occur as intermediate states in vrvr — t- NN, and will 
therefore suppress the isospin index. In these conventions, the Born terms are given by 



9kny 



y=A,s 



Svr {pth)-^ 



Y^imy - m) + my)Qj{y) - m5jQ } , 



with 



_ _t- 2M|, + 2{m\ - m^) 



The partial waves in the partial-wave expansion of the KK scattering amplitude 

oo 

R^'{s,t) = 167r^(2J + l)r^*(t)Pj(cos0f^) 

obey the unitarity relation 



(B.IO) 
(B.ll) 

(B.12) 



j=o 



Imr^*(t)=ar {qthYg'jit) 



JJt 



(t - + CJ, 



K 



+ 



ta 



N 



1 



Hi{t) + Hi{t) 



0[t-tK) 



where gj (t) denotes the vrvr — )• KK partial waves, the KK — )• NN amplitudes fulfill 
and Cj'^ 

I = lb basis (cf. ()B.4p ). With the S'-matrix elements 



H-i{t) = ^{ptktyhi{t), 
pt 



J — 1/2 is an isospin factor that emerges from the conversion between the It 



(B.13) 

(B.14) 
0, 1 and the 



[5j WlifX^ifi. = ^^i^-^/j\t)Q{t-tK) 
we finally obtain the unitarity relation 



^Hi{t)9{t-tN) , (B.15) 



Imhiit) = aj'{r'j{t)yhi{t) 9{t - tK) + <^Iq^' {g^Jit))* fi{t) 9{t - U) , (B.16) 

V2cj 

with TTTT — )■ NN partial waves f±{t) and isospin factors cj for vrA^ scattering (see fT] for precise 
definitions) . 



22 



References 

[1] C. Ditsche, M. Hoferichter, B. Kubis, and U.-G. Meifiner, JHEP 1206 (2012) 043 [arXiv: 1203.4758] 
[hep-ph]]. 

[2] T. P. Cheng and R. F. Dashen, Phys. Rev. Lett. 26 (1971) 594. 

[3] L. S. Brown, W. J. Pardee, and R. D. Peccei, Phys. Rev. D 4 (1971) 2801. 

[4] V. Bernard, N. Kaiser, and U.-G. Meifiner, Phys. Lett. B 389 (1996) 144 |arXiv:hep-ph/9607245] . 

[5] T. Becher and H. Leutwyler, JHEP 0106 (2001) 017 | hep-ph/01032 63l. 

[6] J. Gasser, H. Leutwyler, and M. E. Sainio, Phys. Lett. B 253 (1991) 260. 

[7] J. F. Donoghue, J. Gasser, and H. Leutwyler, Nucl. Phys. B 343 (1990) 341. 

[8] G. Hohler, Pion-Nukleon-Streuung: Methoden und Ergebnisse, in Landolt-Bornstein, 9b2, ed. 
H. Schopper, Springer Verlag, Berlin, 1983. 

[9] J. J. de Swart, M. C. M. Rentmeester, and R. G. E. Timmermans, PiN Newslett. 13 (1997) 96 
[nucl-th/9802084] . 

[10] R. A. Arndt, W. J. Briscoe, I. I. Strakovsky, and R. L. Workman, Phys. Rev. C 74 (2006) 045205 
[nucl-th/0605082] . 

[11] V. Baru, C. Hanhart, M. Hoferichter, B. Kubis, A. Nogga, and D. R. Philhps, Phys. Lett. B 694 
(2011) 473 [arXiv:1003.4444l [nucl-th]]. 

[12] V. Baru, C. Hanhart, M. Hoferichter, B. Kubis, A. Nogga, and D. R. Phillips, Nucl. Phys. A 872 
(2011) 69 [ arXiv:1107.5 509 [nucl-th]]. 

[13] K. Nakamura et al. [Particle Data Group], J. Phys. G 37 (2010) 075021. 

[14] K. M. Watson, Phys. Rev. 95 (1954) 228. 

[15] B. Moussallam, Eur. Phys. J. C 14 (2000) 111 jarXiv:hep-ph/9909292| . 

[16] R. Omnes, Nuovo Cim. 8 (1958) 316. 

[17] S. Descotes-Genon, Thesis, Universite de Paris-Sud, 2000. 

[18] N. I. Muskhelishvili, Singular Integral Equations, Wolters-Noordhoff Publishing, Groningen, 1953 
[Dover Publications, 2nd edition, 2008]. 

[19] P. Biittiker, S . Descotes-Genon, and B. Moussallam, Eur. Phys. J. C 33 (2004) 409 
[aFXiv:hep-ph/031 0283| . 

[20] M. Hoferichter, D. R. Philhps, and C. Schat, Eur. Phys. J. C 71 (2011) 1743 |arXiv:1106.4T47] 
[hep-ph]]. 

[21] I. Caprini, G. Colangelo, and H. Leutwyler, Eur. Phys. J. C 72 (2012) 1860 |arXiv:1111.7T60] 
[hep-ph]]. 

[22] I. Caprini, G. Colangelo, and H. Leutwyler, in preparation. 



23 



[23] D. H. Cohen, D. S. Ayres, R. Diebold, S. L. Kramer, A. J. Pawlicki, and A. B. Wicklund, Phys. 
Rev. D 22 (1980) 2595. 



[24 

[25 
[26 
[27' 
[28 

[29 
[30 
[31 
[32 
[33 
[34; 
[35 

[36 
[37; 

[38 
[39 

[4o; 

[41 
[42; 

[43 
[44; 
[45; 
[46; 
[47; 
[48 
[49; 



A. Etkin et al, Phys. Rev. D 25 (1982) 1786. 

A. V. Anisovich et al, Nucl. Phys. A 662 (2000) 319 far Xiv:1109. 11881 [hep-ex]]. 

R. A. Arndt, W. J. Briscoe, I. I. Strakovsky, and R. L. Workman, Eur. Phys. J. A 35 (2008) 311. 



SAID, jhttp : / / gwdac . phys . gwu . edu 

F. Huang, A. Sibirtsev, J. Haidenbauer, S. Krewald, and U.-G. Meifiner, Eur. Phys. J. A 44 
(2010) 81 [ a?Xiv:0910.4275 ^ [nucl-th]]. 

J. S. Hyslop, R. A. Arndt, L. D. Roper, and R. L. Workman, Phys. Rev. D 46 (1992) 961. 

H. Hellmann, Einfiihrung in die Quantenchemie, Franz Deuticke, Leipzig, 1937. 

R. P. Feynman, Phys. Rev. 56 (1939) 340. 

J. Gasser and H. Leutwyler, Annals Phys. 158 (1984) 142. 

G. Colangelo et al, Eur. Phys. J. C 71 (2011) 1695 [arXiv: 101 1.4408) [heo-lat]]. 
J. Gasser and H. Leutwyler, Nucl. Phys. B 250 (1985) 465. 

B. Ananthanarayan, L Caprini, G. Colangelo, J. Gasser, and H. Leutwyler, Phys. Lett. B 602 
(2004) 218 |arXiv:hep-ph70409222l . 

J. A. Oiler and L. Roca, Phys. Lett. B 651 (2007) 139 j arXiv:0704.0 039] [hep-ph]]. 

R. Garcia-Martm, R. Kamihski, J. R. Pelaez, J. Ruiz de Elvira, and F. J. Yndurain, Phys. Rev. 
D 83 (2011) 074004 [arXiv:1102.2"T83] [hep-ph]]. 

B. Moussallam, Eur. Phys. J. C 71 (2011) 1814 fa rXiv: 11 10.6074 1 [heo-ph]]. 



G. Colangelo, J. Gasser, and H. Leutwyler, Nucl. Phys. B 603 (2001) 125 [arXiv:hep-ph/0103088|. 



M. Frink, B. Kubis, and U.-G. Meifiner, Eur. Phys. J. C 25 (2002) 259 arXiv:hep-ph/0203193| 



J. Bijnens and P. Dhonte, JHEP 0310 (2003) 061 [h ep-ph/0307044| . 

T. Becher and H. Leutwyler, Eur. Phys. J. C 9 (1999) 643 |arXiv:hep-ph/9901384] . 

H. Pagels and W. J. Pardee, Phys. Rev. D 4 (1971) 3335. 

J. Gasser and H. Leutwyler, Phys. Rept. 87 (1982) 77. 

J. Gasser, M. E. Sainio, and A. Svarc, Nucl. Phys. B 307 (1988) 779. 



V. Bernard, N. Kaiser, and U.-G. Meifiner, Z. Phys. C 60 (1993) 111 hep-ph/9303311]. 



R. A. Arndt, R. L. Workman, I. I. Strakovsky, and M. M. Pavan, l arXiv:nucl-th/9807087f 
B. Holzenkamp, K. Holinde, and J. Speth, Nucl. Phys. A 500 (1989) 485. 
W. R. Frazer and J. R. Fulco, Phys. Rev. 117 (1960) 1603. 



24 



